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ABSTRACT 

A model describing Ising spins with short range interactions moving randomly 
in a plane is considered. In the presence of a hard core repulsion, which prevents the 
Ising spins from overlapping, the model is analogous to a dynamically triangulated 
Ising model with spins constrained to move on a flat surface. It is found that as a 
function of coupling strength and hard core repulsion the model exhibits multicritical 
behavior, with first and second order transition lines terminating at a tricritical 
point. The thermal and magnetic exponents computed at the tricritical point are 
consistent with the exact two-matrix model solution of the random Ising model, 
introduced previously to describe the effects of fluctuating geometries. 
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1 Introduction 



Following the exact solution of the Ising model on a random surface by matrix 
model methods [1], there has been a growing interest in the properties of random 
Ising spins coupled to two-dimensional gravity. More recently, work based on both 
series expansions [2] and numerical simulations [3, 4] has verified and extended 
the original results. It is characteristic of these Ising models that the spins are 
allowed to move at random on a discretized version of a fluid surface. In a specific 
implementation of the model, Ising spins are placed at the vertices of a lattice built 
out of equilateral triangles, and the lattice geometry is then allowed to fluctuate 
by varying the local coordination number through a "link flip" operation which 
varies the local connectivity [3]. Remarkably the same critical exponents have also 
been found using consistency conditions derived from conformal field theory for 
central charge c = | [5], which should again apply to Ising spins. It is generally 
believed that the new values for the Ising critical exponents are due to the random 
fluctuations of the surface (or the world sheet in string terminology) in which the 
spins are embedded, and therefore intimately tied to the intrinsic fractal properties 
of fluctuating geometries. It came therefore as a surprise that non-random Ising 
spins, placed on a randomly fluctuating geometry but with fixed spin coordination 
number, exhibited the same critical behavior as in flat space, without any observed 
"gravitational" shift of the exponents [6]. 

The natural question is then to what extent the values of the critical Ising ex- 
ponents found in the matrix model solution (a = —1, (3 = 1/2, 7 = 2, r\ = 2/3, 
v = 3/2 [1]) are due to the annealed randomness of the lattice, and to what extent 
they are due to the physical presence of a fluctuating background metric. The most 
straightforward way to answer this question is to investigate the critical proper- 
ties of annealed random Ising spins, with interactions designed to mimic as closely 
as possible the dynamical triangulation model, but placed in flat two-dimensional 
space. We should add that it is well known that for a quenched random lattice the 
critical exponents are the same as on a regular lattice [7], as expected on the basis 
of universality, even though in two dimensions the Harris criterion (which applies 
to quenched impurities only) does not give a clear prediction, since the specific heat 
exponent vanishes, a = 0, for Onsager's solution. 

In this letter we present some first results concerning the exponents of such a 
model. A more detailed account of our results will be the subject of a forthcoming 
publication [11]. 



2 



2 Formulation of the Model 



In a square d- dimensional box of sides L with periodic boundary conditions we 
introduce a set of N = L d Ising spins Si = ±1 with coordinates x", i = 1...N, 
a = l...d, and average density p = N/L d = 1. Both the spins and the coordinates 
will be considered as dynamical variables in this model. Interactions between the 
spins are determined by 

I[x, S] = - E Mxi, Xj ) Wij SiSj -h^WiSi , (2.1) 

i<j i 

with ferromagnetic coupling 

{0 if | X{ — x j | > R 
J if r < \xi — Xj\ < R , (2-2) 
oo if | X{ — x j | < r 

giving therefore a hard core repulsion radius equal to r/2. As will be discussed 
further below, the hard core repulsive interaction is necessary for obtaining a non- 
trivial phase diagram, and mimics the interaction found in the dynamical triangula- 
tion model, where the minimum distance between any two spins is restricted to be 
one lattice spacing. For r — > 0, Jij = J[l — 8(\xi — xj\ — R)]. The weights Wij and 
Wi appearing in Eq. (2.1) could in principle contain geometric factors associated 
with the random lattice subtended by the points, and involve quantities such as the 
areas of the triangles associated with the vertices, as well as the lengths of the edges 
connecting the sites. In the following we will consider only the simplest case of unit 
weights, Wij = Wi = 1. On the basis of universality of critical behavior one would 
expect that the results should not be too sensitive to such a specific choice, which 
only alters the short distance details of the model. 

The full partition function for coordinates and spins is then written as 

2 = II E (A ft*') exp(-/[z,S]) . (2.3) 

i=l S i = ±l a=l J0 

In the following we will only consider the two-dimensional case, d = 2, for which 
specific predictions are available from the matrix model solution. It should be clear 
that if the interaction range R is of order one, then, for sufficiently large hard 
core repulsion, r — > y/5/2 < R, the spins will tend to lock in into an almost regular 
triangular lattice. As will be shown below, in practice this crossover happens already 
for quite small values of r. The critical behavior is then the one expected for the 
regular Ising model in two dimensions, namely a continuous second order phase 
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transition with the Onsager exponents. Indeed for the Ising model on a triangular 
lattice it is known that J c = |\/31n3 = 0.9514.... On the other hand if the hard 
core repulsion is very small, then for sufficiently low temperatures the spins will tend 
to form tight ordered clusters, in which each spin interacts with a large number of 
neighbors. As will be shown below, this clustering transition is rather sudden and 
strongly first order. Furthermore, where the two transition lines meet inside the 
phase diagram one would expect to find a tricritical point. 

In order to investigate this issue further, we have chosen to study the above 
system by numerical simulation, with both the spins and the coordinates updated 
by a standard Monte Carlo method. The computation of thermodynamic averages 
is quite time consuming in this model, since any spin can in principle interact with 
any other spin as long as they get sufficiently close together. As a consequence, a 
sweep through the lattice requires a number of order N 2 operations, which makes it 
increasingly difficult to study larger and larger lattices. On the other hand, we should 
add that we have not found any anomalous behavior as far as the autocorrelation 
times are concerned, which remain quite comparable to the pure Ising case. 

In the course of the simulation the spontaneous magnetization per spin 

was measured (here the averages involve both the x and S variables, < > = < > Xt s)-> 
as well as the zero field susceptibility 

* = hw^ z ^ = h < ^ SiSi> -h < ^ Si \ >1 ■ (2 - 5) 

It is customary to use the absolute value on the r.h.s., since on a finite lattice the 
spontaneous magnetization, defined without the absolute value, vanishes identically 
even at low temperatures. In addition, in order to determine the latent heat and the 
specific heat exponent, we have computed the average Ising energy per spin defined 
here as 

Id 1 
E = ~N~Qj lnZ \ h=0 = "T/y < J ij( x i, x j) w ij s i s j > > ( 2 - 6 ) 

and its fluctuation, 
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3 Results and Analysis 



In the simulations we have investigated lattice sizes varying from 5 2 = 25 sites 
to 20 2 = 400 sites. The length of our runs varies in the critical region (J ~ J c ) 
between 1M sweeps on the smaller lattices and 100k sweeps on the largest lattices. 
A standard binning procedure then leads to the errors reported in the figures. 

As it stands, the model contains three coupling parameters, the ferromagnetic 
coupling J, the interaction range R and the hard core repulsion parameter r. We 
have fixed R = 1; comparable choices should not change the universality class. As 
we alluded previously, for small r we find that the system undergoes a sharp first 
order transition, between the disordered phase and a phase in which all spins form 
a few very tight magnetized clusters. On the other hand, for sufficiently large r, 
the transition is Ising-like, between ordered and disordered, almost regular, Ising 
lattices (for our choice of range i£, the transition appears to be very close to regular 
Ising-like for r 0.6 and larger, see below). 

A determination of the discontinuity in the average energy of Eq. (2.6) at the 
critical coupling J c shows that it gradually decreases as r is increased from zero. 
Fig. 1 shows a plot of the latent heat versus r at the transition point J c . In general 
we do not expect the latent heat to vanish linearly at the endpoint, but our results 
seem to indicate a behavior quite close to linear. From the data we estimate that 
the latent heat vanishes at r = 0.344(7), thus signaling the presence of a tricritical 
point at the end of the first order transition line. Beyond this point, the transition 
stays second order, as will be discussed further below. The phase transition line 
extends almost vertically through the phase diagram; for r = we found on the 
largest lattices J c = 0.19(2), while for r = 0.6 we found J c = 0.93(3). 

To determine the critical exponents, we resort to a finite size scaling analysis. In 
the following we will be mostly concerned with the values for the critical exponents 
in the vicinity of the tricritical point. In the case of the spin susceptibility, from 
finite-size scaling, we expect a scaling form of the type 

X (N,J) = N^x(N 1/2 »\J- J c \) . (3.1) 

To recover the correct infinite volume result one needs x(x) ~ a; -7 for large argu- 
ments. Thus, in particular the peak in % should scale like jV 7 / 2l/ for sufficiently large 
N. In Fig. 2 we show the evolution of the computed peaks in % as a function of 
IniV, for r = 0.35. 

Despite the fact that the lattices are quite small, as can be seen from the graph 
a linear fit to the data at the tricritical point is rather good, with relatively small 
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deviations from linearity, % 2 /cLo./. = 8.8 X 10~ 3 . Using least-squares one estimates 
7/1/ = 1.27(7), which is much smaller than the exact regular Ising result 7/1/ = 1.75. 
From scaling one then obtains the anomalous dimension exponent r\ = 2 — 7/1/ = 
0.73(7). To further gauge our errors, we have computed the same exponent in the 
regular Ising limit, for r = 0.6. In this case we indeed recover the Onsager value: we 
find on the same size lattices and using the same analysis method 7/1/ = 1.70(8). We 
also note that the shift in the critical point on a finite lattice is determined by the 
correlation length exponent u, namely J C (N) — J c (oo) ~ N~ 1 ^ 2v . This relationship 
can be used to estimate i>, but it is not very accurate. From a fit to the known 
values of J C (N) we obtain the estimate u = 1.3(2). 

A similar finite size scaling analysis can be performed for the magnetization. 
Close to and above J c we expect M ~ ( J — J c )^ , and, at the critical point on a finite 
lattice, as determined from the peak in the susceptibility (which incidentally is very 
close to the inflection point in the magnetization versus J), M should scale to zero 
as Mn(J c ) ~ N^l 2 " ' . In Fig. 3 we show the magnetization M computed in this way 
for different size lattices close to the tricritical point. In spite of the larger errors 
the results again clearly exclude the pure Ising exponents, and give f3/u = 0.30(10), 
to be compared to the exact regular Ising result f3/u = 0.125. A similar analysis in 
the pure Ising limit (more precisely, for r = 0.6) gives f3/u = 0.15(7). 

The results for the peak in the Ising specific heat C at the tricritical point as 
a function of lattice size L are shown in Fig. 4. One expects the peak to grow as 
C ~ jV a / 2l/ , but the absence of any growth implies that a/u < (a weak cusp in 
the specific heat). If we insist on fitting the peak in the specific heat to a power of 
iV, we get a/u —0.11(5), a negative value due to the decrease of the peak with 
increasing system size. On the other hand we should add that, in general close to a 
critical point, the free energy can be decomposed in a regular and a singular part. In 
our case the singular part does not seem to be singular enough to emerge above the 
regular background, leading to an intrinsic uncertainty in the determination of an 
a < 0, and which can only be overcome by determining still higher derivatives of the 
free energy with respect to the coupling J. A better approach would seem therefore 
to determine the correlation length exponent u instead, and use scaling to relate it 
to a = 2 — 2u. In the regular Ising case one has in a finite volume a logarithmic 
divergence C ~ In N (and a/2u = 0), and we indeed see such a divergence clearly 
for r = 0.6, which corresponds to the almost regular triangular Ising case. 

One can improve on the estimate for J c by considering the fourth-order cumulant 
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ll v 


(3/u 


a/u 


V 


This work 


1.27(7) 


0.30(10) 


< -0.11(5) 


1.3(2) 


Matrix model 


1.333... 


0.333... 


-0.666... 


1.5 


Onsager 


1.75 


0.125 


0.0 


1.0 


Tricritical Ising 


1.85 


0.075 


1.60 


0.555... 



Table 1: Estimates of the critical exponents for the random two-dimensional Isin^ 
model, as obtained from finite size scaling at the tricritical point. 
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where m 



U N {J) 



< > 



(3.2) 



3 < m 2 > 2 ' 

S{/N. It has the scaling form expected of a dimensionless quantity 



U N (J) = U(N 1 ^\J- J c \). 



(3.3) 



The curves Un(J), for different and sufficiently large values of iV, should then in- 
tersect at a common point J c , where the theory exhibits scale invariance, and U 
takes on the fixed point value U* . We have found that indeed the curves meet very 
close to a common point, and from the intersection of the curves for N = 25 to 
400 we estimate J c = 0.48(1), which is consistent with the estimate of the critical 
point derived from the location of the peak in the magnetic susceptibility. We also 
determine U* = 0.47(4), to be compared to the pure Ising model estimate for the 
invariant charge U* 0.613 [10]. 

One can estimate the correlation length exponent u from the scaling of the slope 
of the cumulant at J c . For two lattice sizes N,N' one computes the estimator 

\n[N'/N] 



"eff(N,N') 



(3.4) 



2\n[U' N ,(J c )/U' N (J c )} 

with U' N = 8Un /dJ. Using this method, we find u = 1.3(3). 

In Table I we summarize our results, together with the exponents obtained for the 
two-matrix model [1], for the Onsager solution of the square lattice Ising model, and 
for the tricritical Ising model in two dimensions [8]. As can be seen, the exponents 
are quite close to the matrix model values (the pure Ising exponents seem to be 
excluded by several standard deviations). 



4 Conclusions 

In the previous sections we have presented some first results for a random Ising 
model in flat two-dimensional space. The model reproduces some of the features 
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of a model for dynamically triangulated Ising spins, and in particular its random 
nature, but does not incorporate any effects due to curvature. Due to the non- 
local nature of the interactions of the spins, only relatively small systems could be 
considered so far, which is reflected in the still rather large uncertainties associated 
with the exponents. Still a rich phase diagram has emerged, with a tricritical point 
separating first from second order transition lines. We have localized the tricritical 
point at J c = 0.48(1) and r = 0.344(7). The thermal and magnetic exponents 
determined in the vicinity of the tricritical point (presented in Table I) have been 
found to be consistent, within errors, with the matrix model solution of the random 
Ising model. Our results would therefore suggest that matrix model solutions can 
also be used to describe a class of annealed random systems in flat space. 
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Figure 1: Latent heat along the first order transition line, plotted against the hard 
core repulsion parameter r. The tricritical point is located where the latent heat 
vanishes. 



Figure 2: Peak in the magnetic susceptibility, Xmaxi versus the number of Ising spins 
iV, for fixed hard core repulsion parameter r = 0.35. 



Figure 3: Finite size scaling of the magnetization at the inflection point, Mi n f, versus 
the total number of Ising spins iV, for fixed hard core repulsion parameter r = 0.35. 



Figure 4: Plot of specific heat C versus ferromagnetic coupling J at r=0.35, showing 
the absence of a growth in the peak with increasing lattice sizes, in contrast to the 
behavior of the magnetic susceptibility. The errors (not shown) are smaller than the 
size of the symbols. 
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